Phase transitions in Ising magnetic films and superlattices 
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lii; Abstract 
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QQ . Within the framework of mean field theory, we examine the phase tran- 
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cn . sitions in Ising magnetic films and superlattices. By transfer matrix 

O 

-(— > 



method, we derive two general nonlinear equations for phase transition 
temperatures of Ising magnetic films and superlattices, respectively. The 
equations can be applied to the films and superlattices with arbitrary 
exchange interaction constants and arbitrary layer number. Numeri- 



o . 

Q , cal results for phase transition temperatures as a function of exchange 

- ^^ , interaction constants are presented. 
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I. INTRODUCTION 

Considerable effort has been recently devoted to the understanding of magnetic 
filnis, layered structures and superlattices[l-7]. With the development of molecular 
beam epitaxy, it is now possible to grow in a very controlled way magnetic films with 
few atomic layers or even monolayer atop nonmagnetic substrates. A superlattice in 
which the atoms vary from one monolayer to another can also be envisaged. Very 
often one finds unexpected and interesting properties in these systems. For exam- 
ple, experimental studies[8-10] on the magnetic properties of surfaces of Gd,Cr and 
Tb have shown that a magnetically ordered surface can coexist with a magnetically 
disordered bulk phase. 

Most work has been devoted to the free surface problem in which spins on the 
surface interact with each other through an exchange interaction J a different from that 
of bulk material [11-20]. For the values of Jg/ J above a certain critical value {Js/J)crit 
the system orders on the surface before it orders in the bulk. Below this critical 
value only two phases are expected, namely the bulk ferromagnetic and paramagnetic 
phases. 

Magnetic excitations in superlattices were studied in numerous papers (see, e.g., 
Ref.21 for a brief review). Yet less attention has been paid to critical behavoir. 
The phase transition temperatures for a Heisenberg magnetic superlattices have been 
studied [22-24]. Hinchey and Mills have investigated a superlattice structure with 
alternating ferromagnetic and antiferromagnetic layers [25-26]. 

In this paper, we are concerned with phase transitions in Ising magnetic films 
and superlattices. we study them within the framework of mean field theory. The 
transfer matrix mathod is used to derive nonlinear equations for magnetic Ising films 
and superlattices. The equations are general ones for arbitrary exchange interaction 
constants. In section II, we outline the formalism and derive the nonlinear equation 
for phase transition temperatures of magnetic films. The transition temperatures 
as a function of surface exchange constants are studied. The nonlinear equation for 



transition temperatures of superlattices are given in section III. The last section IV 
is devoted to a brief discussion. 

II. FORMALISM AND PHASE TRANSITIONS IN ISING MAGNETIC 

FILMS 

We start with a lattice of localized spins with spin equal to 1/2. The interaction 
is of the nearest-neighbor ferromagnetic Ising type. The Ising Hamiltonian of the 
system is given by 

H = —- 2_^ 2_^ JijSirSjr', (1) 

{i,j) {r,r') 

where {i,j) are plane indices, {r,r') are different sites of the planes, and Sir is spin 
variable. Jij denote the exchange constants and are plane dependent. We will keep 
only nearest-neighbour terms. 

In the mean-field theory. Sir is replaced by its mean value rrii associated with each 
plane, and is determined by a set of simultaneous equations 

rrii = tanh[{zoJiimi + zJi^i+inii^i + zJi^i_imi_i)/KBT], (2) 

where Zq, z are the numbers of nearest neighbours in the place and between the planes, 
respectively, /c^ is the Boltzman constant and T is the temperature. 

Near the transition temperature, the order parameter rrii are small, Eq.(2) reduces 
to 

KBTrrii = ZoJanii + zJj^j+imj+i + zJi.i_imi_i. (3) 

Let us rewrite the above equation in matrix form in analogy with Ref. [27] 

> =M.-i('"-) (4) 

\ rrii J \mi-ij 

with Mj„i as the transfer matrix defined by 

(\KbT — ZQJii)/[zJi^i^i) —Ji,i-ii Ji,i+i \ 
(5) 



We consider a magnetic film which contains A^ layers with layer indices i = 1, 2...N. 
From Eq.(4), we get 

f'"M=flh), (6) 

where R = Mjv-2---^2^i represents successive multiplication of the transfer matrices 

For an ideal film system, there exists symmetry in the direction perpendicular to 
the surface, which allows us to write rrii = -mTv+i-j- Then, the following nonlinear 
equation for the transition temperature can be obtained from equations (3) and (6) 
as 

i?n [{KbT - zoJii)/izJ^,2)? + {Ri2 - R2i)[{KbT - z,J^^/{zJ,^^] - R22 = 0.(7) 

The above equation is the general equation for the transition temperature of sym- 
metric films. It is suitable for arbitrary exchange interaction constants Jij. All the 
information about the phase transition temperatures of the system is contained in 
the equation. 

For a uniform system with J^ = J, Eq.(7) reduces to 

Rn{t - zo/zf + {R12 - i?2i)(t - zo/z) - R22 = 0, (8) 

where t = KbT/{Jz) is the reduced temperature and R = D^"'^. Here the matrix 
ft-zo/z -1\ 

Note that det(-D) = 1, we can linearalize matrix R as[28] 

R = Un^2A - U^.^I, (10) 

where / is the unit matrix ,f/jv = (A^ — A^)/(A+ — A_), and \± = (t — zq/z ± 



(t — Zq/z)^ — 4)/2. Substituting Eq.(lO) into Eq.(8), we reduce Eq.(8) to its sim- 
plest form 

f/iv+i = 0. (11) 



f/jv+i can be rewritten as 

f/iv+i = sin[(iV+l)0]/sin0 (12) 

for (t — zq/ zY < 4. Here = arccos[(t — z^/ z)/2]. For (t — zq/ zY > 4, becomes 19, 
and the trigonometric functions become hyperbolic functions of 9. 
Eq.(ll) gives 

t = 2cos[7r/(A^ + l)]+ V-2- (13) 

In the hmit A^ — > oo, the reduced bulk temperature ts is obtained as 

tB = 2 + z^/z. (14) 

Fig.l shows the reduced transition temperature vs. layer number L. It can be seen 
that the film transition tempeature is lower than the bulk one, i.e., the film disorders 
at a lower temperature than the bulk ones. Throughout this paper, we take J as the 
unit of energy. 

In the above discussions, we have assumed the surface exchange constants Js are 
the same as the bulk exchange constants J . Now we consider a (/, n, /) film consisting 
/ top surface layers, n bulk surface layers and / bottom surface layers, and assume that 
the exchange constants in a surface layer is denoted by Js and that in a bulk layer or 
between successive layers by J . In this case the total transfer matrix R becomes 

R = P^~^Q''P^~\ (15) 

where the matrix P and Q are 

(t-AJs/J -l\ (t-A ~l\ 

Here we have assumed that the spins lie on a simple cubic lattice, i.e., Zq = 4, z = 1. 
Since Det{P) = Det{Q) = 1, the matrix p'^^ and Q"' can be linearalized in analogous 
to Eq.(lO). In this case, the nonlinear Eq.(7) reduces to 

i?n(t - AJs/Jf + (i?i2 - i?2i)(t - AJs/J) - R22 = 0. (17) 



The numerical results for transition temperatures of a magnetic {l,n,l) film as a 
function of Js/J are shown in Fig.2. It is can be seen that the transition tempera- 
ture increase as Js/J increases. For Js/J < 1, the transition temperatures in film 
(10, 10, 10) are nearly equal to the bulk temperature Tb- It is interesting that the 
transtion temperature increases linearly with the increase of Js/J when Jg/J is large 
enough. We can also see that the transition temperature increases as layer number 
increases. 

III. PHASE TRANSITIONS IN ISING MAGNETIC SUPERLATTICES 

The (/,n) superlattice structure we study is formed from two types of atoms. In 
each elementary unit with layer indices i = 1,2. ..I + n, there are / atomic layers of 
type A and n atomic layers of type B. The interlayer exchange constants are given 
by Ja and Jf,, whereas the exchange constants between different layers is described by 
J. For the above model of the superlattice, the transfer matrix Mj(Eq.(5)) reduce to 
different types of matrix 

-(? ;')-(? ■; 

where Xa = t—AJA/J and Xb = t—AJB/J- From Eq.(4), we can obtain the following 
equation as 

h---A=RM, (19) 

where 

R = AB'^A^-^ (20) 

is the total transfer matrix. 

Due to the periodicity of the superlattice, we know mi^n+2 = 1^2 and rrii^n+i = 
mi. Then from Eq.(19), we get 

Det{R) - Tr{R) + 1=0. (21) 



It can be easily seen that Det{A) = Det{B) = Det{R) = 1. Then the Eq.(21) reduces 
to the simplest form 

Tr{R) = 2. (22) 

Actually, the above equation is a general equation for phase transition temperature 
of superlattices. It is valid for arbitrary exchange constants Jij. The nonlinear eqution 
for films are dependent on both diagonal and nondiagonal terms of the total transfer 
matrix R, while the nonlinear equation for superlattices are only depend on diagonal 
terms of R. 

For the total transfer matrix R = AB^''A''^^, we get 

TriA^B"") = 2. (23) 

The matrix A'' and -B" can be linearlized as 

A' = EiA - Ei^J 
E" = F^B - Fn^J, (24) 

where E^ = «-«'_)/(«+-«-), F„ = {f3l-f3^)/{f3^-f3_), a± = (X^±y^Xl - 4)/2 



and P± = {Xb ± yX^ — 4)/2. From Eq.(23), we get the equation 

(EiXa — Ei^i){FnXi) — F„„i) — 2EiFn + Ei„iFn-i = 2. (25) 

For / = 1, n = 1, the above equation reduces to 

XaX, = 4, (26) 

which is identical with the result of Ref.(2) and (29). Next we numerically calculate 
the phase transition temperatures from Eq.(25). 

In figure 3, we have shown the results for the superlattices (3, 1), (2, 2), (10, 5) 
and (20,20). The transition temperature is plotted as a function of Ja/J- For 
J a/ J < 1, the transition temperature is smaller that the bulk transition temper- 
ature. For Ja/J = 1, the transition temperature of the superlattice is independent of 
m and n. and are equal to the bulk temperature Tg as expected. On the other hand. 



for Ja/ J > Ij the transition temperature is greater than the bulk temperature Tb- 
The transition temperature increases with the layer number in one unit cell and ap- 
proaches Tb asymptotically as the number become large. The transition temperature 
increases nearly linearly with Ja/J when Ja/J is large enough. The layer number of 
superlattices (3, 1) and (2,2) are same, but the transition temperatures are different. 
For Ja/J < 1, the transition temperatures of superlattice (2, 2) are larger than those 
of superlattice (3, 1). In contrary to this, the transition temperatures of superlattice 
(2, 2) are smaller than those of superlattice (3, 1) for Ja/J > 1- 

IV. DISCUSSIONS 

In summary, we have studied phase transitions in Ising magnetic films and su- 
perlattices with the framework of mean field theory. By transfer matrix method, we 
have derived two general nonlinear equations for phase transition temperatures in 
Ising films and superlattices, respectively. The transition temperatures as a function 
of exchange interaction constants are calculated. In addition, the equations can be 
easily solved and the parameters involved can be adjusted at will. 



Figure Captions 

Fig.l, Transition temperatures of a uniform film as a function of layer number L. 
Fig. 2, Transition temperatures of a magnetic film (/, n, /) as a function of Js/J- 
Fig. 3, Transition temperatures of a magnetic superlattice (/, n) as a function of Ja/J- 
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